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Abstract 

As yet, certain aspects of the Africanization process are not well understood, for example, the reproductive behavior 
of African and European honeybees and how the first Africanized swarms were formed and spread. Drone congrega- 
tion areas (DCAs) are the ideal place to study honeybee reproduction under natural conditions since hundreds of 
drones from various colonies gather together in the same geographical area for mating. In the present study, we as- 
sessed the genetic structure of seven drone congregations and four commercial European-derived and Africanized 
apiaries in southern Brazil, employing seven microsatellite loci for this purpose. We also estimated the number of 
mother-colonies that drones of a specific DCA originated from. Pairwise comparison failed to reveal any population 
sub-structuring among the DCAs, thus indicating low mutual genetic differentiation. We also observed high genetic 
similarity between colonies of commercial apiaries and DCAs, besides a slight contribution from a European-derived 
apiary to a DCA formed nearby. Africanized DCAs seem to have a somewhat different genetic structure when com- 
pared to the European. 
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Introduction 

Since the introduction of the African honeybee sub- 
species Apis mellifera scutellata into Brazil in 1956 there 
has been a dramatic change in the genetic makeup of Euro- 
pean-derived honeybee populations in South and Central 
America in less than 50 years (Schneider et al, 2004). The 
replacement of the European by Africanized populations, 
and the strong predominance of African traits in the hybrid 
fornied with the spread of African swarms throughout the 
New World, have been well documented (Lobo et al. , 1 989; 
Smith etal., 1989; Quezada-Euan and Paxton, 1999; Clarke 
et al. , 200 1 ; Collet et al , 2006). 

Although several reasons for the prevalence of Afri- 
can traits in the resultant hybrid have been proposed 
(Schneider et al., 2004), two important issues are, as yet, 
not well understood, namely, the reproductive behavior of 
African and European honeybees and how the first Afri- 
canized swarms were formed and spread. The reproductive 
behavior of ^. mellifera under natural conditions can be as- 
sessed through studies on drone congregation areas 
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(DCAs). These are where hundreds of drones from various 
colonies gather in the same geographical area for mating. 
The main properties of a DCA are: 1) the highly mixed 
drone population from many surrounding apiaries; 2) the 
formation of these areas takes place independent of the 
presence of queens; and 3) the clustering of flying drones 
can be found every year in the same place (Jean-Prost, 
1958). The composition of drone congregations is of great 
consequence for the genetic structure of honeybee colonies, 
as the number of colonies represented in a DCA affects the 
genetic diversity of new colonies (Baudry et al, 1998). 
Thus, DCAs are ideal for studying honeybee reproduction 
under natural conditions. 

Most of the studies concerning DCA properties were 
hereto based on honeybee populations in Europe (Baudry et 
al, 1998; Koeniger et al, 2005a, 2005b; Kraus et al. 
2005a), Asia (Punchihewa et al, 1990; Koeniger et al, 
1994; Kraus et al, 2005b; Wattanachaiyingcharoen et al, 
2008) or Africa (Free, 1987; Muerrle et al, 2007). Few 
studies on Africanized DCAs have been undertaken, these 
being related to drone racial identification for European 
queen-mating inferences in African-derived honeybee ar- 
eas (Loper and Fierro, 1991), to the identification of the de- 
gree of Africanization in certain areas (Eischen and 
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Rubink, 1997), or yet to assess the seasonal abundance of 
European and African-derived drones in DCAs for areas 
undergoing Africanization (Quezada-Euan and May-ltza, 
2001). There are no studies concerning the genetic structure 
of DCAs in Africanized areas of South and Central Amer- 
ica. 

In southern Brazil, near the towns of Lajeado and 
Eldorado do Sul (Rio Grande do Sul State), seven DCAs, 
located close to commercial apiaries of European-derived 
or Africanized colonies, were identified. The presence of 
European-derived and Africanized colonies in the area pre- 
sented an opportunity for checking how these bees cooper- 
ate in forming a DCA. In this study, and by means of seven 
microsatelite loci, we analyzed the genetic structure of 
seven drone congregations and four commercial Euro- 
pean-derived and Africanized apiaries. Besides estimating 
the number of mother colonies drones originated from, we 
also analyzed whether those visiting the DCAs originated 
from one or more subpopulations. Herein, we will describe 
(1) the low genetic differentiation among seven DCAs, (2) 
the genetic similarity between each DCA and the surround- 
ing apiaries, and (3) the low representation of the European 
mitotype in those DCAs close to European-derived colo- 
nies. 

Material and Methods 
Sampling 

A total of 410 ^. melUfera drones were sampled from 
seven DCAs located near Africanized (T4, T6, T7 and SI) 
or European-derived apiaries (CM2, CM3, and HC) in 
Lajeado (30° 05' 44.1" S; 51° 40' 07.8" W) and Eldorado 
do Sul (29° 25' 43.1" S; 51° 58' 29.5" W) (Figure 1). All 
the DCAs were sampled in September, 2002, when the pro- 
duction of drones is high. Samples from CM2 were col- 
lected on two subsequent days and at different times, and 
were denominated CM2H1 (14:00-15:00 h) and CM2H2 
(15:00-16:00 h). All drone congregations were very close, 
one to the other, except for DCA HC, which was located 
76 km apart (Figure 1). Flying drones were captured with 
an aerial trap (Taylor, 1984) and kept at -20 °C until use. 
Workers were collected from four apiaries (Funcionarios, 
Tambo, Casa-mel and Hans Cremer) located near these 
DCAs, so as to assess the effects of the presence of Afri- 
canized and European-derived colonies on DCA composi- 
tion. The numbers of drones and workers per DCA and 
apiary are given in Table 1 . 

DNA extraction and analysis 

Total DNA was extracted from the legs of the drones 
or from the thorax of a worker (one worker per colony) 
from each apiary by means of the phenol-chloroform 
method (Sheppard and McPheron, 1991). 

Seven microsatellite loci (A14, A24, A43, A79, A88, 
A113 e B124) were amplified according to Estoup et al. 
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Figure 1 - Location of the DCAs and apiaries sites in Lajeado and Eldo- 
rado do Sul (Rio Grande do Sul state). 



(1995) and Solignac et al (2003). Polymerase chain reac- 
tion (PCR) amplifications were performed in a 1 0 |iL total 
volume containing 1 U of Platinum Taq DNA polymerase 
(Invitrogen), lOx Taq DNA polymerase buffer, 2.2 mM 
MgCb (except for A24, A43 and B124 with 2.5 mM 
MgClz), 100 |iM of each dNTP, 0.5 |aM of each primer and 
1 |iL of template DNA for all primer pairs. The forward 
primer for each marker was labeled with one of three fluo- 
rescent dyes (Fan, Vic or Ned; Applied Biosystems). PCR 
conditions consisted of an initial 2 min denaturation step at 
94 °C followed by 30 cycles of 94 °C for 30 s, annealing 
temperatures at 55 °C (A24, A43 and B 124) or 62 °C (A 14, 
A79, A88 e Al 13) for 30 s, and 72 °C for 30 s, plus a final 
extension at 72 °C for 10 min. The PCR products were run 
on a MegaBace 750-GE capillary sequencer using ET 
550-R MegaBace (GE) as the internal size standard. 
MegaBace Fragment Profiler 1 .2 software was used for al- 
lele identification. 

The mtDNA pattern of the samples from the apiaries 
and DCAs was identified through Dral restriction patterns 
of the COI-COII intergenic region for the former, and Dral 
restriction patterns of the 16S fragment for the latter, re- 
spectively. The polymerase chain reaction (PCR) was car- 
ried out according to Gamery et al. (1992) for COI-COII, 
and digestion with Dral enzyme to Collet et al. (2006). 
Amplification and restriction of the 16S fragment was un- 
dertaken according to Collet et al. (2007). 

Data analysis 

The number of alleles per locus (n) and expected 
heterozygozity (He) were estimated with GENEPOP v3.4 
software (Raymond and Rousset, 1995). We used either 



Genetics of Africanized drone congregation areas 



859 



Fstat v2.9.3 software (Goudet, 2001) to calculate pairwise 
Fst values and their significant deviation from zero, as well 
as allelic richness (a^), or the sum of the number of alleles 
actually observed in each population when considering dif- 
ferences in sample size. In order to reveal whether there is 
genetic structure encompassing the data set, molecular 
variance analysis (AMOVA) was undertaken with Arle- 
quin v2000 software (Schneider et al, 2000). The same 
program was used to assess hierarchical multi-loci 
AMOVA. Variance was then partitioned into variance 



components and distributed among the groups (HC + T4, 
T6, T7, SI + CM2, CM3), and among and within DCAs. 

Colony vl .3 software (Wang, 2004) was used for esti- 
mating the number of colonies contributing drones to a 
given DCA. Five replicate analyses were performed with 
different seed numbers to calculate the number of colonies. 

Results 

The number of alleles, allelic richness and genetic di- 
versity of the DCAs and apiaries appear in Table 1. The 
number of alleles ranged from 5 (A24) to 17 (A14), with an 



Table 1 - Sample size in each location (N), number of alleles (n), allelic richness (a, ) and genetic diversity (Hd) by microsatellite loci in Apis mellifera 
drone congregation areas and apiaries from southern Brazil. 
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average of 11.15 + 2.2 alleles per locus for DCAs and 
11.03 + 2.4 for the apiaries. We found an average allelic 
richness of 7.9 +1.3 alleles per locus for the DCAs and 
8.3 + 1.5 for the apiaries. 

Each drone sampled had its mtDNA origin assessed 
through Dra\ restriction patterns of the 16S fragment. The 
DCA T7 and HC showed one and two drones, respectively, 
with a European pattern (carnica/ligustica) while all the re- 
maining drones (n = 407) exhibited an African pattern. 
Dml restriction patterns of the COI-COII mtDNA region 
permitted classifying the samples from apiaries located 
near the DCAs as either European (E), Africanized (A) or 
mixed (E/A) (Table 2). 

No significant difference (Fst = 0.0025) was found 
between samples of the DCA CM2 collected during two 
consecutive days in two different periods. Due to the ge- 

Table 2 - Classification of the samples from apiaries located near the 
DCAs as European (E), Africanized (A) and both (E/A) according to 
COI-COII (Dral) patterns. 
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netic homogeneity, CM2 was considered as a single unit for 
further analyses. 

The DCA spatial pattern was reflected at the genetic 
level. Hierarchical AMOVA revealed that the greatest vari- 
ation (99.9%) occurred among individuals within the 
DCAs themselves (Table 3). As a consequence, the six 
DCAs at Eldorado do Sul, closely distributed along ~5 km, 
did not show any population sub-structuring. The pairwise 
comparison also showed a low differentiation among most 
of the DCAs and apiaries with European-derived colonies 
(Figure 1 and Table 4). As expected, Hans Cremer and 
Casa-mel apiaries showed the greatest differentiation to the 
DCAs sampled (Table 4). 

Based on genotypic data, we estimated the number 
of colonies contributing drones to each DCA as; T4- 
17.8 ± 0.45; T6-26.2 ± 0.84; T7-26; Sl-8; CM2-58 ± 1.23; 
CM3-13; HC - 21.2 ± 0.45. Since the F,t values revealed 
low levels of differentiation, we considered all DCAs as 
constituting a single unit, in order to estimate the number 
of drone-contributing colonies. The total number of colo- 
nies was 100.6 ± 1.34 (mean ± SD), thus proving that 
DCAs are constituted of poorly related drones of 
multi-colonial origin. The distribution of drones over the 
colonies followed a normal distribution (Wilks Shapiro 
Normality Test, W = 0.91; p > 0.1), ranging from one to 
eight contributed drones per colony, without any sign of 
skewed distribution or over-representation from any sin- 
gle colony (Figure 2). 



Table 3 - Analysis of Molecular Variance (AMOVA) conducted for drone congregation areas of Apis mellifera. (df = degree of freedom; SS = sum of 
squares). 

Source of Variation df SS Variance components % variation p-value 

Among groups 2 0.933 -0.00049 -0.10 > 0.05 

Among DCAs within groups 6 3.124 0.00055 0.11 > 0.05 

Within DCAs 401 199.621 0.49781 99.99 > 0.05 



Table 4 - Pairwise Fji values between pairs of Apis mellifera drone congregation areas. 





T4 


16 


T7 


SI 


CM3 


HC 


CM2 H. Cremer 


Func. 


Tambo 


T4 




















T6 


-0.0030 


















T7 


-0.0055 


-0.0011 
















SI 


0.0007 


-0.0112 


-0.0084 














CM3 


-0.0059 


-0.0034 


0.0018 


-0.0044 












HC 


-0.0009 


0.0021 


-0.0025 


-0.0082 


0.0022 










CM2 


0.0064 


0.0072 


0.0079* 


-0.0155 


-0.0015 


0.0040 








H. Cremer 


0.0462* 


0.0561* 


0.0457* 


0.0479 


0.0488* 


0.0310 


0.0414* 






Func. 


-0.0018 


0.0033 


-0.0026 


-0.0022 


0.0004 


-0.0057 


-0.0027 0.0352* 






Tambo 


0.0019 


-0.0013 


0.0008 


-0.0028 


0.0038 


-0.0107 


-0.0142 0.0361 


0.0040 




Casa-mel 


0.0408* 


0.0301* 


0.0281* 


0.0175 


0.0292 


0.0192 


0.0319* 0.0594* 


0.0227 


0.0194 



*p<0.05. 



Genetics of Africanized drone congregation areas 



861 



40 




I 2 3 4 5 6 7 8 

Number of drones 

Figure 2 - Number of colonies contributing drones to the DCAs. 



Discussion 

Drones from the Africanized DCAs showed a higher 
number of alleles and gene diversity in comparison to a 
DCA from Gemiany (A. m. carnica) studied by Baudry et 
al. (1998) (Table 5). Higher gene variation, expressed as 
the expected heterozygosity, was found in the Africanized 
DCAs. African populations are significantly more variable 
than European, this having been interpreted as a conse- 
quence of the higher density of colonies and higher migra- 
tory ability of African subspecies when compared to the 
Europeans (Estoup etal, 1995; Francke^a/., 1998). More- 
over, Africanized honeybees are hybrid populations. Thus, 
a higher degree of genetic variation due to the contribution 
of some European and African alleles to the gene pool of 
the hybrid Africanized bee would be expected. 

The most frequent alleles found in European popula- 
tions (Estoup et al. 1995, Franck et al. 1998) are also the 
most frequent among Brazilian DCAs. This can be inter- 
preted as either a significant European contribution to Afri- 
canized populations, or that these alleles are not efficient 
racial markers, since they are shared by the subspecies in- 
volved in the Africanization process. The latter interpreta- 
tion seems to be more feasible, as allozymic (Lobo et al, 
1989; Del Lama et al, 1990), mtDNA (Smith et al, 1989; 
Clarke et al, 2001; Collet et al, 2006) and morphometric 



data (Diniz-Filho and Malaspina, 1995) all indicate a very 
high proportion of A. m. scutellata genes in Africanized 
honeybees. 

No significant differences were found between sam- 
ples of the DCA CM2 collected during two different peri- 
ods. This finding is in accordance with previous observa- 
tions that honeybees from European-derived or Africanized 
colonies do not present flight-activity at different periods in 
Neotropical regions (Cristino AS, unpublished results). 

In Brazil, we found a higher variation in the number 
of drones contributed by a colony to a DCA, following a 
normal distribution. In contrast, most colonies in European 
DCAs (75%) contribute with the same number of drones, 
following a Poisson distribution (Baudry et al, 1998). This 
differential contribution in Africanized colonies could be 
explained by differences in colony drone production due to 
genetic differences. The different ability for drone- 
production in the first African swarms would result in 
strong selection pressure, and therefore must have in- 
creased the differential drone contribution of African and 
European colonies to DCAs. Certainly, this was not the sole 
factor to determine the observed degree of mixture that 
tends towards African mitotypes. To address this question 
properly, it would be necessary to account for the additive 
genetic variability linked to drone production in Africa- 
nized colonies, this being an important trait to be worked by 
artificial selection. 

Three out of the 410 drones analyzed showed the 
carnicalligustica mtDNA pattern. The under-represen- 
tation of European drones in a DCA was particularly evi- 
dent in DCA HC, located 500 m from the Hans Cremer 
apiary (Figure 1 and Table 3), with colonies aggregating a 
large number of drones. This apiary was initially consti- 
tuted of 40 carnica queens brought from Europe, the mater- 
nal lineage having been maintained through the production 
of new carnica queens. This result did not confirm the 
statement that drones choose the nearest DCA and that this 
choice boosts the genetic representation of local colonies 
(Koeniger et al, 2005a). 



Table 5 - Comparison between Brazilian and European samples for the number of alleles (nj and expected heterozygozity {H^. For the Brazilian sam- 
ples, we present the average number of alleles and H^. 
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The distribution of maternal nuclear alleles in each 
DCA and its closest apiary was not significantly different. 
These observations are consistent only if we consider that 
males copulating with a queen are a representative sample 
of the local bee population. The high number of colonies 
represented in a drone congregation makes mating between 
a queen and one of her brothers unlikely. If many colonies 
contribute equally to the DCA, a queen has an equal chance 
of insemination by a male from each colony within the re- 
cruitment area of a DCA. This results in low average relat- 
edness between a queen's mates, thereby avoiding high 
inbreeding and homozigosis at the csd locus (Beye et al, 
2003; Beye, 2004) for her progeny, and consequently pre- 
venting the production of diploid drones, so avoiding a pos- 
sible, heavy genetic load for the colony. Thus, high 
polyandry associated with the composition of drone con- 
gregations ensure very high intra-colony gene variation, 
thereby resulting in an efficient panmitic system. 

The estimated number of colonies contributing with 
drones to the DCAs was much lower, when compared to the 
number presented by Baudry et al. (1998) for a carnica 
DCA. This result can be explained by the lower number of 
microsatellite loci used in this study, thereby decreasing the 
chances of detecting the number of families of brothers in 
each DCA. However, the most likely explanation for this 
result must be the distinct process by which DCAs are 
formed in Europe/USA and Brazil. Whereas DCAs are well 
defined and in small number in Europe and the USA, in 
Brazil they seem to be more continuously distributed. At 
least eight DCAs were identified in a study carried out in an 
area with a radius of three kilometers in Ribeirao Preto (Sao 
Paulo State, Brazil) (Martinez-Caranton OA and Scares 
AEE, unpublished results). Moreover, while in Europe the 
colonies are concentrated in apiaries and only a few feral 
swarms are found, in Brazil the number of feral swarms is 
high. Therefore, if a higher number of areas exist, one could 
expect a lower number of drones per DCA in Brazil and, as 
a consequence, a lower number of contributing colonies. 

The major results of this study are the low genetic dif- 
ferentiation among the DCAs, the genetic similarity be- 
tween workers of commercial apiaries and DCAs fonned 
nearby, and the small contribution of a European apiary to a 
DCA fornied nearby. The lack of genetic differentiation in 
the studied DCAs was also observed in samples of apiaries 
from 43 localities from northern to southern Brazil for 
twelve microsatellite loci (unpublished data). The homoge- 
neous genetic composition of the DCAs showed on a 
micro-scale what happens on a macro-scale, i.e., the homo- 
geneous genetic conditions of Brazilian Africanized bee 
populations. Considering that 26 scutellata queens colo- 
nized much of the Americas and largely replaced European 
bees throughout their range in the New World, it is very im- 
pressive how this founder effect produced what has been 
considered one of the prime examples of a successful bio- 
logical invasion known. 
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